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Some time ago a novel method of investigation of fluctuations in even-by-event 
analysis of high energy multiparticle production data has been proposed ^. It 
consists in defining a suitable measure $ of a given observable x being exactly 
the same for nucleon-nucleon and nucleus-nucleus collisions if the later are simple 
superpositions of the former. 
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where Z = '^ Zi. (1) 



i=l 



Here Zi = Xi — x where x denotes the mean value of the observable x calculated for 
all particles from all events (the so called inclusive mean) and A^ is the number of 
particles analysed in the event. In (|I|) (A^) and {Z"^) are averages of event-by-event 
observables over all events whereas the last term is the square root of the second 
moment of the inclusive z distribution. By construction |I| if particles are produced 
independently ^^^ = 0. 

When applied to the recent NA49 data from central Pb — Pb collisions at 158 
A-GeV this method revealed that fluctuations of transverse momentum (x = px) 
decreased significantly in respect to elementary NN collisions. This has been inter- 
preted as a possible sign of equilibration taking place in heavy ion collisions provid- 
ing thus enviroment for the possible creation of quark-gluon plasma (QGP). It was 
quickly realised that existing models of multiparticle production are leading in that 
matter to conflicting statements [0. On the other hand the use of fluctuations as a 
very sensitive tool for the analysis of dynamics of multiparticle reactions has been 
advocated since already some time [Q , especially their role in searching for some spe- 
cial features of the QGP equation of state has been shown to be of special interest |^. 

However, it was demonstrated recently in that the corresponding fluctuation 
measure calculated for a pion gas in global equilibrium (defined within the standard 
extensive thermodynamic) is almost an order of magnitude greater then the exper- 
imental value. Although the recent NA49 paper presents a new value, which 
is more like the prediction in 0, the controversy aroused around $ resulted in a 
number of presentations trying to clarify and extend the meaning of the $ variable 
(cf., for example, [§] and references thereinQ)^. In the mean time the use of this 

-•^For some other recent discussions of event-by-event fluctuations see [^. 

^We would like to point here only that, if there are some additional fluctuations (not arising 
from quantum statistics, like those caused by the experimental errors) which add in the same way 
to both terms in definition (|I|) of $ it would perhaps be better to use $ ^- $* = -^^^ — z^ -where 
they would cancel. 



variable has been extended, so far only theoretically, to the possible study (actu- 
ally planned already by NA49) of the event-by-event fluctuations of the "chemical" 
(particle type) composition of the final stage of high energy collisions [0, |Tl[ . 



Moving apart from the above discussion we would like to follow another path 



of research. Namely, it was suggested recently in |T^ that the extreme conditions 
of density and temperature occuring in ultrarelativistic heavy ion collisions can 
lead to memory effects and long-range colour interactions and to the presence of 
non-Markovian processes in the corresponding kinetic equations (cf., for example 
[[13| ). It turns out that such effects in many other branches of physics are best 
described phenomenologically in terms of a single parameter q by using the so called 
nonextensive statistics |]n| . This statistics is based on new definition of g-entropy 



(which for g — > 1 coincides with the usual Boltzmann-Gibbs definition): 

^ w w 

Sq = —— 51 Pfe (l - pV^) '^^ S = -Y. Pk InPfc (2) 

^ ^ fc=i fe=i 

(defined for the probability distribution {p^} for a system of W microstates). Such 

entropy is nonextensive, i.e., for a system {A + B) composed with two independent 

systems A and B: 

S,{A + B) = S,iA) + S,{B) + (1 - g) Sg{A) ■ S,{B). (3) 

The value of |g — 1| characterises deviation from the extensitivity, see [12, l3 for 



more details. As was shown in |T^, nonextensive approach with q as minuite as 
q = 1.01 -T- 1.015 eliminates the abovementioned discrepancy between the first NA49 
data and ideal quantum gas (Boltzmann statistics) estimation of 0. 

Because NA49 Collaboration plans to study also the chemical fluctuations, there 
have already apeared predictions concerning the expected form of the fluctuation 
measure $ in this case |jlll. They are based on the use of normal (i.e., Boltzmann) 
statistics. In the present note we have generalized it to the case of nonextensive 



statistics in a manner identical to that presented in |12| for the case of transverse 



momenta. As in |TT| we have computed the $ measure for the system of particles of 
two sorts, 7r~ and K~ , i.e., nonstrange and strange hadrons with multiplicities (ri^) 
and {uk), respectively. Since 

{N) = {n^) + (uk) (4) 

one immediately finds that in definition (|l]) 

^' ' ~ (AT) ^^^ 



and 

(^ ) - ^vp • ^^> 

He have now consistently replaced the mean occupation numbers by their g-equivalents, 
which under some approximations, valid for small values of nonextensitivity |1 — g|, 



can be expressed in the following analytical form |15 



(n), = {[1 + (g - imE - ^t^"^'-^^ ± l}"' , (7) 

where /3 = 1/kT, /i is chemical potential and the +/— sign applies to fermions/bosons. 
Notice that in the limit g ^ 1 (extensive statistics) one recovers the conventional 
Fermi-Dirac and Bose-Einstein distributions. We shall not dwell on the details of 



this procedure, they are essentially the same as those discussed in [0. It is only 
necessary to mention that in this approximation one retains the basic factorised 
formula for correlations used in |TT[, namely that 



{niTij) = {ni){nj). (8) 

As in |1^ (where $ for transverse momenta pt has been considered), a rather large 



sensitivity of predictions presented in [O to the parameter q has been observed. 



According to the nonextensive statistics philosophy this fact indicates a large sensi- 
tivity to the (initial and boundary) conditions present in the ultrarelativistic heavy 
ion collisions and existence of some kind of memory effects in such systems, as 
mentioned in references [0. Our results are presented in Figs. 1 and 2 where 



modifications caused by the nonextensitivity q = 1.015 (chosen in such a way as to 
fit the Pt spectra in Fig. 3, see discussion below) to the results of [0 for directly 



produced particles are shown. For simplicity we have restricted ourselves here only 
to comparison with results of |jll| without resonances R. 



Notice (cf. also [0) that, as is clearly seen in Fig. 3, the same pattern of 



fluctuations is already present in the transverse momentum spectra of produced sec- 
ondaries, i.e., the same value of q brings new "g-thermal" curve in agreement with 
experiment in the whole range of pt presented. If it would emerge also in the future 
data on the fluctuations of chemical composition discussed here, i.e., if parameter 
q would turn out to be similar (modulo experimental errors), it would signal that 

^One can argue that resonance production belongs in our philosophy already to the nonextensive 
case being therefore responsible for (at least a part of) the effect leading to a nonzero |1 — g]. This 
is best seen inspecting results of [Q with resonances included, which show that $ in this case also 
changes sign. The use of parameter q is, however, more general as it includes all other possible 
effects as well. 



both observables, fluctuations of wliicli is investigated, are similarly affected by the 
external conditions mentioned before and that they can be easily parametrized phe- 
nomenologically by a single parameter g, i.e., by the measure of the nonextensitivity 
of the nuclear collision process^. 



^It is worth to mention here that methods of nonextensive statistics have been already used in 
the field of high energy physics in order to analyse some aspects of cosmic ray data fl^ and to the 



description of hadronization in e"'"e~ annihilation processes |18 . Also the recently proposed use of 
quantum groups in studying Bose-Einstein correlations observed in all multiparticle reactions ||l^ 
belong to that cathegory because, as was shown in pfl], there is close correspondence between the 
deformation parameter of quantum groups and the nonextensitivity parameter of Tsallis statistics. 
In fact, as can be seen from |1J], also works on intermittency using the so called Levy stable 
distributions (for example pTl) belong to this cathegory as well. 
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Figure Captions: 



Fig. 1 $ - measure of the kaon multiplicity fluctuations (in the n^K^ system of 
particles) as a function of temperature for three values of the pion chemical 
potential. The kaon chemical potential vanishes. The resonances are neglected, 
(a) - results of |Tl[ (in linear scale); (6) - our results for q = 1.015. 



Fig. 2 $ - measure of the kaon multiplicity fluctuations (in the 7r~if~ system of 
particles) as a function of temperature for three values of the kaon chemical 
potential. The pion chemical potential vanishes. The resonances are neglected, 
(a) - results of |Tl| (in linear scale); (6) - our results for q = 1.015. 



Fig. 3 The results for px distribution: notice that q = 1.015 results describes also the 
tail of distribution not fltted by the conventional exponent (i.e., g = 1 in our 



case, cf. also ||T2|)- Data are taken from |jT6| 
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